function y = func_s_from_p(p, CONSTS)

    eps_a = CONSTS.eps_a;
    y = (p.^2 - eps_a).^0.5; 
    
end
